quant-ph /040 1 1 54 



Optimal Database Search: Waves and Catalysis 

Apoorva Patel 

Centre for High Energy Physics, Indian Institute of Science, Bangalore-560012, Indi^ 

Grover's database search algorithm, ahhough discovered in the context of quantum computation, 
can be implemented using any system that allows superposition of states. A physical realization of 
this algorithm is described using coupled simple harmonic oscillators, which can be exactly solved in 
both classical and quantum domains. Classical wave algorithms are far more stable against decoher- 
ence compared to their quantum counterparts. In addition to providing convenient demonstration 
models, they may have a role in practical situations, such as catalysis. 
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I. THE OPTIMAL SEARCH ALGORITHM 

Database search is an elementary computational task 
with wide-ranging applications. Its efficiency is measured 
in terms of the number of queries one has to make to 
the database in order to find the desired item. In the 
conventional formulation of the problem, the query is a 
binary oracle (i.e. a YES/NO question). For an unsorted 
database of N items, using classical Boolean logic, one 
requires on the average {Q) = N binary queries to locate 
the desired item. The number of queries is reduced to 
(Q) = (A^ + l)/2, if the search process has a memory 
so that an item rejected once is not picked up again for 
inspection. 

Grover discovered a search algorithm that, using super- 
position of states, reduces the number of required queries 
to Q = O^sfN) pj]. This algorithm starts with a super- 
position state, where each item has an equal probability 
to get picked, and evolves it to a target state where only 
the desired item can get picked. Following Dirac's nota- 
tion, and using the index i to label the items, the starting 
and target state satisfy 
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The algorithm evolves \s) towards by discrete rota- 
tions in the two-dimensional space formed by \s) and 
using the two reflection operators. 



Ut^l-2\t){t\ , C/, = l-2|s)(s| 
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Ut is the binary oracle which flips the sign of target 
state amplitude, while —U^ performs the reflection-in- 
the-mean operation. Solution to Eq.(3) determines the 
number of queries as 



(4) 



(2g + l)sin-^(l/VAf) = 7r/2 



(In practice, Q must be an integer, while Eq.(4) may not 
have an integer solution. In such cases, the algorithm is 
stopped when the state has evolved sufflciently close to, 
although not exactly equal to, \t). Then one finds the 
desired item with a high probability.) 

In the qubit implementation of the algorithm, one 
chooses N = 2^ and the items in the database are la- 
beled with binary digits. Using a uniform superposition 
as the starting state. 
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(i|s) = l/VA^, [/, ==H«"(l-2|0)(0|)ff^ 
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[H is the Hadamard operator), the implementation re- 
quires only 0(log2 A^) spatial resources [1]. It has been 
proved that this is the optimal algorithm for unsorted 
database search 

Several variations of this algorithm have been studied, 
using different physical realizations of the database items 
and the target query oracle. In the original algorithm, 
the states are encoded in an rt-qubit register, and the ora- 
cle is a discrete binary operation (denoted by Ut above). 
In the analogue version of the algorithm, the discrete 
unitary oracle is traded for a continuous time interaction 
Hamiltonian (it evolves the target state somewhat dif- 
ferently than the rest) that acts for the entire duration 
of the algorithm, and the number of queries is replaced 
by the time one has to wait for before finding the target 
state Q. The wave version of the algorithm requires A^ 
distinct wave modes, instead of n qubits, but does not 
involve quantum entanglement at any stage Such a 
wave search has also been experimentally implemented 
using classical Fourier optics, with a phase-shift plate 
providing the oracle Q. A classical coupled pendulum 
model of the analogue version of the algorithm has been 
described, where one of the pendulums is slightly differ- 
ent from the rest and the uniform superposition state \s) 
is identified with the center-of-mass mode 0]. In what 
follows, I describe a binary oracle version of the wave 
search algorithm using identical coupled harmonic oscil- 
lators. 
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II. HARMONIC OSCILLATOR 
IMPLEMENTATION 

A system of coupled harmonic oscillators is frequently 
studied in physics. It involves only quadratic forms, and 
hence can be solved exactly in both classical and quan- 
tum domains. This property makes it extremely useful in 
situations where cross-over between classical and quan- 
tum behaviour is to be analyzed. We shall first look at 
the classical system, and then observe that the quantum 
system essentially follows the same pattern. 



A. Classical oscillators 

Let the items in the database be represented by N 
identical harmonic oscillators. While they are oscillating 
in a specific manner, someone taps one of the oscillators 
(i.e. elastically reflects it by touching it). The task is to 
identify which of the oscillators has been tapped, without 
looking at the tapping. The optimization criterion is to 
design the system of oscillators, and their initial state, so 
as to make the identification as quickly as possible. 

Grover's algorithm requires identical coupling between 
any pair of oscillators. This can be accomplished by cou- 
pling all the oscillators to a big oscillator, as shown in 
Fig.l. The big oscillator then becomes an intermediary 
between any pair of oscillators, with the same strength, 
since it is coupled to the center-of-mass mode. The La- 
grangian for the whole system is 
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£ = -MX^ - -KX^ + Y.[^mx^ - ^Hx, - Xf] . (6) 



With the center-of-mass displacement, x = Y^f^iXi/N , 
the Lagrangian can be rewritten as 
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Now we can fix the oscillator parameters to implement 
Grover's algorithm. In the algorithm, we are interested 
in the dynamics of the tapped oscillator. All the other 
oscillators (i.e. i ^ t) influence the dynamics of xt only 
through the combination x. The dynamics of [N — 2) 
linearly independent modes orthogonal to Xt and x (they 
all have the form {xj^t —Xk^t)) decouples from the modes 
of interest; we can drop them and effectively work in the 
3-dimensional space of the modes xt}. (In what 

follows, we shall first specify initial conditions such that 
all Xi^t are identical and all the decoupled modes vanish. 
Subsequently, we will look at the general situation by 
adding back all the decoupled modes.) 

Ghoosing units of mass and time such that m = 1, fc = 
1, and in terms of the variables 

Y^VmX, y^VNx, yt^xt-x, (8) 

the effective Lagrangian becomes 
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The potential energy terms in £off are easily diagonalized, 
and yield the eigenvalues 




The corresponding eigenmodes are 

e± = (1 - Lul)Y + = (1 - loDVmX + -^x , 
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The initial uniform superposition state can be realized 
as all the oscillators moving together, while the big oscil- 
lator is at rest. 
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t = 0: X ^0 , X ^0 , x,^Q 
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(We will consider situations with general initial condi- 
tions later.) The reflection operators correspond to shift- 
ing the appropriate oscillator phases by half a period. 



[a) 

(b) 
(c) 



Xt 
Xt 



0, Xt > : Just before the oracle 
0, Xt —Xt : Elastic reflection 
0, if < : Just after the oracle 



FIG. 1: A system of identical harmonic oscillators, coupled 
to a big oscillator via the center-of-mass mode. 



FIG. 2: The binary tapping oracle flips the sign of the target 
oscillator velocity, when its displacement is zero. 
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The binary tapping oracle can be realized as the elastic 
reflection illustrated in Fig. 2. That implements Ut in the 
velocity space, by reversing the target oscillator velocity 
at the instance when all the displacements vanish. Time 
evolution of the coupled oscillators redistributes the total 
kinetic energy, and that can implement the operator Us 
with a suitable choice of time interval and frequencies. 

Wit h the natural frequency of individual oscillators 
uj — ^/kjrn — 1, the reflection- in-the- mean operation re- 
quires uj± to be rational numbers. Optimization means 
that they should be selected to make the dynamics of the 
whole system of oscillators have as small a period as pos- 
sible. The solution is not unique. One set of solutions is 
(p is a positive integer) 



M = 



2p+l 



16Nm 



3{2p + 3)i2p ~ 1) 



K 



3{2p + 3){2p-i 



In these cases, in the absence of oracles, the dynamics 
of the whole system of oscillators has the period, T = 
Att. The big oscillator returns to its initial rest state 
{X = 0,X = 0), whenever t is an integral multiple of 
27r, i.e. after every half a period. Time evolution for the 
same interval of half a period reverses a;, while leaving 
it — X unchanged, i.e. it implements the operator Us in 
the velocity space. Thus Grover's algorithm, Eq.(3), can 
be realized by applying the tapping oracle at every time 
interval At = 27t. 

A more interesting set of solutions is 

Nth 

a;+ = 2p , w_ = =^ M = , if = . 

+ ^' (2p+l)(2p-l) ' 

(14) 

In these cases, the big oscillator is not coupled to any sup- 
port, and e_ becomes a translation mode for the whole 
system of oscillators. The translation mode can be elim- 
inated from the dynamics with the initial conditions 

N 

t=:0: X = , X = -— v4 , = , = A . (15) 

Then, in the absence of oracles, the dynamics of the 
whole system of oscillators has the smallest possible pe- 
riod, T = 2tt. After half a period, the big oscillator is 
back to its initial state, x also returns to its initial value, 
while it — X changes its sign. This is equivalent to ap- 
plying —Us in the velocity space, and Grover's algorithm 
can be implemented by tapping the target oscillator at 
every time interval At — tt. 

There is an important physical distinction between the 
quantum and the wave interpretations of the amplitude 
amplification process in Grover's algorithm — quantum 
probability is mapped to wave energy. The enhancement 
of the quantum amplitude increases the probability of 
finding the target state iV-fold, while the enhancement 
of the wave amplitude increases the energy of the tar- 
get oscillator A^-fold. The well-known phenomenon of 



"beats" is responsible for energy transfer amongst cou- 
pled oscillators. The elastic reflection oracle does not 
change energy, and it is interesting to observe that the 
oscillator which is obstructed by tapping picks up energy. 



B. Stability considerations 

Now we can look at the behaviour of the wave imple- 
mentation under general circumstances. First consider 
the initial conditions. Despite appearances, precise syn- 
chronization of oscillators is not an issue in the algorithm, 
because of the explicit coupling to the center-of-mass 
mode. For instance, the algorithm can be started off 
with an initial push to the big oscillator, X = B,ii = 0, 
and the system of oscillators would evolve to the stage 
X = 0, Xi = A. Furthermore, any arbitrary distribution 
of initial velocities of oscillators can be accommodated 
in the analysis by bringing back the {N — 2) decoupled 
modes. The decoupled modes have no effect whatsoever 
on the dynamics of the {X,x,xt} modes. Consequently, 
the algorithm is only modified to the extent that the en- 
ergy amplification of the target oscillator is limited to the 
initial energy present in the {X^x,Xt} modes, instead of 
being iV-fold. Explicitly, the maximum gain is 



Nx 



N 



N- 1 



{it 



ix? 



(16) 



which can be substantial for the generic situation where 
the initial it and x are of the same order of magnitude. 

To extract the maximum gain, the algorithm must be 
stopped at a precise instant (i.e at a precise value of 
Q); otherwise the evolution continues in repetitive cycles. 
The state evolution in Grover's algorithm is a uniform 
rotation in the two dimensional |s)-|t) subspace. The av- 
erage overlap of the target state, with the state \q{Q)) 
after Q queries, is therefore 



l(<z(g)|i)lL = (sin'^?)av = l/2. 



(17) 



Thus if the algorithm is stopped at a random instant, the 
energy gain on the average is half of its maximum value 
in Ea. ((T6|) — which can still be substantially larger than 
1. 

Next consider the effect of damping The crucial 
ingredient in the algorithm is the coherence amongst the 
phases of the oscillators. That is governed by the fre- 
quencies of the oscillators, and is independent of the am- 
plitudes. For a weakly damped oscillator, its amplitude 
changes linearly with the damping coefficient, while its 
frequency changes quadratically. The time evolution of 
the above implementation, therefore, remains essentially 
unaffected if the oscillators experience a small damping. 
The leading effect is a decrease in the energy amplifica- 
tion due to decaying amplitudes. 

Among other variations, simultaneous scaling of 
masses and spring constants of the oscillators (i.e. rrii = 
aim and ki = Uik) does not alter the algorithm at all. 
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since the scale factors can be absorbed by redefining Xi. 
One can also contemplate some changes of global condi- 
tions that let the algorithm go through but change its 
physical interpretation. For example: (i) Application of 
simultaneous tapping oracle to more than one oscillator 
(with suitable changes in M and K) can focus the energy 
into the tapped oscillators, (ii) Interchange of initial and 
final states can run the algorithm backwards, whence the 
large initial energy of the target oscillator gets uniformly 
distributed among all oscillators, (iii) Elastic reflection 
of all but the target oscillator implements the oracle —Ut 
in the velocity space, in which case it is the unobstructed 
oscillator that picks up energy. 



C. Quantum domain 

The dynamics of harmonic oscillators is simple enough 
to permit exact quantum analysis as well. It is convenient 
to interpolate between classical and quantum domains 
using the coherent state formulation fsl . Coherent states 
are superpositions of the eigenstates, parametrized by a 
single complex variable a. 



|a) 
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of the half-oscillator are the same as those for the har- 
monic oscillator, with odd n. It is straightforward to 
ensure the node at x = using the method of images, 
and the tapped oscillator coherent states become 
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(24) 



with the normalization constant C — (1 — e 



-2|q|^-j-1/2 

Tapping amounts to change-over between |a) and \ — a), 
which reverses (x) and (p) compared to the untapped mo- 
tion. In addition, the wavefunction changes sign, which 
is the geometric phase corresponding to wave reflection. 
The evolution of a coherent state wavepacket undergoing 
reflection from the wall is depicted in Fig. 3. 



III. POSSIBLE APPLICATIONS 

The oscillator based search process discussed above 
has the same algorithmic efficiency as the proposals of 
Refs.jl, H, 0] — where it differs from them is in the ac- 
tual physical implementation. Wave algorithms are clas- 
sical, but they have not been explored much historically. 
The main reason is that they require exponentially more 
spatial resources compared to their digital counterparts. 



They describe Gaussian wavepackets with minimal 
spread (i.e. displaced versions of the ground state eigen- 
function). 
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A coherent state with the initial condition a{t = 0) = 
has energy Tiwdaop + j), and the centre of its wavepacket 
performs the same simple harmonic motion as a classical 
particle would: 



{^)(t) , (p){t) 



(20) 
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The wavefunction of the state evolves according to 

|V(0)) = |ao) =^ |7/^(i)) =e— ^ (21) 
while the explicit structure of the wavepacket is given by 

1/4 /x~(x)\^ .x{p) 
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(22) 



Thus the analysis of Section II. A can be carried over 
unchanged to the quantum domain, provided we can fig- 
ure out how the tapping oracle works for the coherent 
states. The tapped oscillator corresponds to a particle 
moving in the half-oscillator potential 

V{x) ^ \kx^ for X < , V{x) = oo for x > . (23) 

The impenetrable wall at x = is equivalent to enforcing 
the boundary condition -(/'(x = 0) = 0. So the eigenstates 
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FIG. 3: Evolution of the coherent state wavefunction of the 
tapped oscillator, with the initial condition qq = —a, The 
left half of the figure shows the actual wavepacket in the 
harmonic oscillator potential, with the impenetrable wall at 
a; = 0. The right half of the figure shows the image wavepacket 
that ensures the node of the wavefunction at x = 0. For 
t — nLj/2, the wavefunction is purely imaginary, but the fac- 
tor of i is omitted for convenience in drawing the figure. The 
wavepacket at t = ttoj includes the geometric phase of — 1 
arising from reflection. 
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0{N) vs. 0(log2 N). On the other hand, they can reduce 
the number of oracle calls by exploiting superposition of 
states. Note that no algorithm based on Boolean logic, 
either with serial or with parallel implementation, can re- 
duce the number of oracle calls for an unsorted database 
search to less than 0{N). 

Quantum algorithms are superior to wave algorithms, 
because they can use superposition as well as reduce spa- 
tial resources. The reduction of spatial resources, how- 
ever, comes with the cost that quantum algorithms have 
to work with entangled states. Quantum entanglement is 
far more sensitive to decoherence caused by environmen- 
tal disturbances than mere superposition, and that has 
made physical implementations of quantum algorithms 
very difficult. On the other hand, superposition of classi- 
cal waves can be fairly stable, even when a small amount 
of damping is present, and that can make wave imple- 
mentations advantageous in specific physical contexts. 

These comparisons suggest that wave algorithms fall in 
a regime inbetween classical and quantum algorithms — 
more efficient than the former and more robust than the 
latter. They are likely to be useful in practical situa- 
tions, where N is not very large and environmental dis- 
turbances are not negligible. Indeed it is worthwhile to 
systematically explore them, just like randomized algo- 
rithms have been 

In the specific case of the unsorted database search 
problem, the remarkable simplicity of the oscillator im- 
plementation makes one wonder about possible appli- 
cations, besides constructing convenient demonstration 
models. The implementation of Rcf. 6] has one oscillator 
with a frequency different from the rest, which gets sin- 
gled out by dynamical evolution of the coupled system. 
In the present scenario, all oscillators are identical, but 
one of them is discretely tapped by an external agency 
and gets selected. I point out a physical situation below, 
that can fit such a scenario, and where involvement of 
new mechanisms can enhance our understanding of the 
observed phenomena. 



A. Catalysis 

The practically useful property of the wave search al- 
gorithm is that it focuses energy into one of the oscillator 
modes. There exist a large number of chemical reactions 
which, although not forbidden by energy conservation, 
are extremely slow because they have to pass through an 
intermediate state of high energy. In these reactions, the 
dominant term governing the reaction rate is the Boltz- 
mann factor, exp(— i?b/fcr), with the barrier energy in 
the exponent. Only a tiny fraction of the molecules in 
the tail of the thermal distribution arc energetic enough 
to go over the barrier and complete the reaction. It is 
known that the rates of many such reactions can be en- 
hanced by orders of magnitude by adding suitable cat- 
alysts (enzymes in case of biochemical reactions) to the 
reactants. The conventional explanation for the reaction 



rate enhancement, called transition state theory, is that 
the catalysts lower the energy of the intermediate state 
by modifying the chemical environment of the reactants. 

The preceding analysis of the wave search algorithm 
suggests another mechanism for catalysis. Vibrations 
and rotations of molecules are ubiquitous harmonic os- 
cillator modes. The catalyst can act as the big oscilla- 
tor and focus energy of many modes into the reactant 
which faces the energy barrier. For example, the cat- 
alytic substrate can have many identical molecules of one 
reactant stuck to it and vibrating, the second reactant 
then comes drifting along and interacts with one of the 
stuck molecules, that molecule picks up energy from its 
neighbours and the reaction gets completed. In such a 
scenario, for maximum efficiency, the physical parame- 
ters (masses and spring constants) need to have specific 
values. But even without perfectly tuned parameters, 
there can be partial energy focusing that provides useful 
increase in the reaction rate. Whether this mechanism 
exists among the known catalysts, or whether we can de- 
sign new type of catalysts that use this mechanism, is an 
open question. 

The catalytic role of chemical environment vs. physi- 
cal waves can be tested by isotopic substitution in the 
reactants, since isotopic substitution changes physical 
parameters without altering chemical properties. The 
conventional transition state theory has no mass depen- 
dence, so any isotope dependence of reaction rates is a 
signal of involvement of physical (in contrast to chemi- 
cal) features in the process. Many examples of isotopic 
dependence of catalytic reaction rates have been discov- 
ered, and the effect is referred to as the "Westheimer ef- 
fect" or the "kinetic isotope effect" [l0[ . The effect is the 
largest for substitution of hydrogen by deuterium, and 
has been extensively studied for the rupture of C-H/C-D 
bonds. The transition state theory has been found inad- 
equate for theoretical understanding of the observations, 
and vibrationally enhanced quantum tunneling has been 
invoked as an alternative [ll|. 

In this context, the oscillator based search process de- 
scribed above has two novel features to suggest yet an- 
other alternative: 

(i) The energy focusing mechanism enhances reaction 
rates beyond their naive classical values. This enhance- 
ment is both mass and temperature dependent, because 
the energy enhancement depends on the masses involved 
while the Boltzmann factor contains the temperature. 

(ii) The wavefunction sign-flip caused by reflection can 
switch between bonding and anti-bonding molecular or- 
bitals, and thus help in the transfer of atoms. This fea- 
ture is not related to either mass or temperature. 

It would be therefore worthwhile to explore the oscilla- 
tor inspired catalytic mechanism with careful modeling 
of specific reactions. 

In closing, I mention an intriguing possibility. I had 
constructed a quantum database search scenario for the 
DNA replication process, whose stumbling block was the 
maintenance of quantum coherence in presence of con- 
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tinuous jostling of molecules [I4I. It should be feasible 
to construct a wave database search scenario involving 
vibrational and rotational modes of the molecules. That 



would be far less troubled by decoherence effects, and 
make interpretation of the genetic language as optimal 
solution to the database search process likely. 
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